Reaction dynamics through kinetic transition states 
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The transformation of a system from one state to another is often mediated by a bottleneck in the 
system's phase space. In chemistry these bottlenecks are known as transition states through which 
the system has to pass in order to evolve from reactants to products. The chemical reactions are 
usually associated with configurational changes where the reactants and products states correspond, 
e.g., to two different isomers or the undissociated and dissociated state of a molecule or cluster. In 
this letter we report on a new type of bottleneck which mediates kinetic rather than configura- 
tional changes. The phase space structures associated with such kinetic transition states and their 
dynamical implications are discussed for the rotational vibrational motion of a triatomic molecule. 
An outline of more general related phase space structures with important dynamical implications is 
given. 
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Introduction. - Reaction type dynamics is characterized 
by the property that a system spends long times in one 
phase space region (the region of the 'reactants') and oc- 
casionally finds its way through a bottleneck to another 
phase space region (the region of the 'products'). The 
main examples are chemical reactions where the bottle- 
necks are induced by saddle points of the potential en- 
ergy surface which arise from a Born-Oppenheimer ap- 
proximation and determine the interactions between the 
constituent atoms or molecules involved in the reaction. 
The reactions are then characterized by configurational 
changes like, e.g., in an isomerization or disscociation 
reaction. In this case the bottleneck is referred to a tran- 
sition state. The main idea of transition state theory 
which is the most frequently used approach to compute 
reaction rates is then to define a surface in the tran- 
sition state region and compute the rate from the flux 
through this so called dividing surface. For getting the 
exact rate this way, it is crucial to define the dividing 
surface in such a way that it is crossed exactly once by 
reactive trajectories (trajectories that evolve from reac- 
tants to products) and not crossed at all by nonreactive 
trajectories (i.e. trajectories which stay in the reactants 
or products region). In the 1970's it has been shown by 
Pechukas, Pollak and others that for systems with two 
degrees of freedom, such a dividing surface can be con- 
structed from an unstable periodic orbit which gives rise 
to the so called periodic orbit dividing surface (PODS) 
[HE]. It took several decades to understand how this 
idea can be generalized to systems with three or more 
degrees of freedom pj]. The object which replaces the 
periodic orbit is a so called normally hyperbolic invariant 
manifold (NHIM) g]. The NHIM does not only allow for 
the construction of a dividing surface with the desired 
crossing properties but it also has (like the unstable peri- 
odic orbit in the two-degree-of-freedom case) stable and 
unstable manifolds which have sufficient dimensionality 



to form separatrices that channel the reactive trajectories 
from reactants to products and separate them from the 
nonreactive ones [5]. In this paper we show that for sys- 
tems which are not of type 'kinetic-plus-potential', there 
are saddle type equilibrium points which either instead 
of the bottlenecks associated with configuration changes 
induce bottlenecks of a kinetic nature or, more gener- 
ally, do not even induce bottlenecks but still give rise to 
NHIMs whose stable and unstable manifolds govern the 
dynamics near such equilibrium points. 

The standard case. - For a system of type 'kinetic-plus- 
potential' (so called natural mechanical systems 0), sad- 
dle points of the potential lead to equilibrium points of 
Hamilton's equations near which the Hamiltonian can be 
brought through a suitable choice of canonical coordi- 
nates (so called normal form coordinates 0[7]) into the 
form 
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where / = n + 1 is the number of degrees of freedom 
and Eq is the energy of the saddle of the potential. 
The quadratic part of the Hamiltonian ([I]) consists of 
a parabolic barrier in the first degree of freedom whose 
steepness is characterized by A > (the Lyapunov expo- 
nent) and n harmonic oscillators with frequencies u>k > 0, 
k = 1, . . . , ri. 

Let us ignore the higher order terms for a moment and 
rewrite the energy equation H = E in the form 

^ + Ef^+^= B - £ o+^o- (2) 
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Then one sees that for E > Eq (i.e. for energies above the 
energy of the saddle), each fixed value of the reaction co- 
ordinate qo defines a (2/ — 2)-dimensional sphere S 2 ?~ 2 . 
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The energy surface thus has the topology of a 'spherical 
cylinder', S 2 ^~ 2 x HL This cylinder has a wide-narrow- 
wide geometry where the spheres S 2 f~ 2 are 'smallest' 
when q — (see [5] for a more precise statement). In fact 
the (2/ — 2)-dimensional sphere given by setting go = 
on the energy surface H = E > Eq defines a dividing sur- 
face which separates the energy surface into a reactants 
region go < and a products region go > 0. All forward 
reactive trajectories cross it with p > 0. All backward 
reactive trajectories cross it with p a < 0. The condition 
Po = defines a (2/ — 3)-dimensional sphere which di- 
vides the dividing surface into the two hemispheres which 
have po > and po < and hence can be viewed as 
the 'equator' of the diving surface. In fact the equator 
is a normally hyperbolic invariant manifold (NHIM) [3]. 
It can be identified with the transition state: a kind of 
unstable 'super molecule' sitting between reactants and 
products [9] . The NHIM has stable and unstable man- 
ifolds S 2f - 2 x R. They thus have one dimension less 
than the energy surface H — E. They form the phase 
space conduits for reaction [10] . The NHIM is of cen- 
tral importance because it dominates the dynamics in 
the neighborhood of the saddle. An important aspect of 
the theory of NHIMs is that they persist if the higher 
order terms in ([TJ) are taken into account and the energy 
is not too far above E [3]. The phase space structures 
above persist accordingly. 

Kinetic transition states. - Let us now more generally 
consider the case of an equilibrium point of a Hamiltonian 
system which has one pair of real eigenvalues ±A and n 
complex conjugate pairs of imaginary eigenvalues iiw/t, 
k = l,...,n. Choosing normal form coordinates near 
such a so called saddle-center-- • - -center equilibrium the 
Hamiltonian assumes the form (JlJ. In fact this general 
case already formed the starting point in [3l [5] and later 
works. However all these studies were restricted to the 
case where the frequencies uik are positive. In this paper 



we study what happens if we give up this assumption. 
To this end let us for simplicity consider a system with 
/ = 2 degrees of freedom with Hamiltonian function 

H = E + \(pl-a 2 ) + U -(p\ + q \) (3) 

with A > and 10 < 0. In his case let us rewrite the 
energy equation H = E in the form 

- 2 <ll-^{pl + <Z) = - 2 pl + E«-E. (4) 

This is the same as ^ however with the role of go and po 

exchanged and E — E Q replaced by E — E. Accordingly 

for E < Eq, the right hand side is positive for any po and 

the energy surface has again the structure of a spherical 

cylinder with a wide-narrow-wide geometry. The crucial 

difference to the standard case of positive frequencies in 

Eq. ([2]) is that the reaction coordinate is now Pq instead of 

go. The bottleneck is thus associated with kinetic rather 
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FIG. 1: Definition of Jacobi coordinates for a triatomic molecule. 

than configurational changes. As the following example 
shows such kinetic bottlenecks do indeed exist in many 
important applications. 

Example: rotational vibrational motion of triatomic 
molecules.- The Hamiltonian of a triatomic molecule is 
given by [HI HI] 
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+ V(pi,P2,#. (5) 



Here (pi,p2,<fi) are Jacobi coordinates defined as 

/0X = ||si||, P2 = |MI) Si • s 2 = pip 2 cos0, 

where si and s 2 are the mass-weighted Jacobi vectors 
(see Fig. [TJ 

, — / \ , — / mixi -r-m 3 x 3 

s l = VAM X 1 - x 3j i S 2 = V7*2( X 2 ; ), 
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computed from the position vectors x& of the atoms and 



the reduced masses 
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The momenta p%, p 2 and p^ in §5§ are conjugate to p\, 
P2 and 4>, respectively, and J = (Ji, J2, J3) is the body 
angular momentum. The magnitude J of J is conserved 
under the dynamics generated by the Hamiltonian ([5]). 

Let us at first consider a rigid molecule (i.e. the val- 
ues of pi, P2 and 4> are fixed). The body fixed frame 



FIG. 3: Potential V as a function of the Jacobi coordinate <p. 



FIG. 2: Angular momentum sphere with contours of the Hamilto- 
nian in j8|. 

can then be chosen such that the moment of inertial ten- 
sor becomes diagonal with the principal moments inertia 
Mi < M 2 < M3 ordered by magnitude on the diagonal. 
The Hamiltonian (fsl) then reduces to 




As the magnitude J of J is conserved the angular mo- 
mentum sphere J 2 + Jf + Jf = J 2 can be viewed as 
the phase space of the rigid molecule [T2] . The solution 
curves of this one-degree-of-freedom system are obtained 
from the level sets of the Hamiltonian H on the angu- 
lar momentum sphere (see Fig. [2]) . The Hamiltonian H 
has local minima at (Ji, J2, J3) = (0, 0, ± J) of energy 
J 2 /(2M 3 ), local maxima at (J u J 2 , J 3 ) = (±J,0,0) of 
energy J 2 /(2Mi) and saddles at (Ji, J 2 , J 3 ) = (0, ± J, 0) 
of energy J 2 / (2M 2 ). These correspond to the center equi- 
libria of stable rotations about the first and third prin- 
cipal axes and the saddle equilibria of unstable rotations 
about the second principal axis, respectively. A possi- 
ble choice of canonical coordinates (q,p) on the angular 
momentum sphere is |12j 

Ji = \/ J 2 - p 2 cos q, J 2 = \/ J 2 -p 2 sing, J 3 =p. (9) 

Since (Ji,J 2 ,Jz) = (J, 0,0) resp. (q,p) = (0,0) is a 
maximum of the Hamiltonian the normal form of the 
Hamiltonian at this equilibrium is H = J 2 /(2Mi) + 
ui(Pq + g 2 )/2+h.o.t. with the negative frequency u> — 
- y/M 2 M 3 (M 2 - Mi)(M 3 - Mi)/MiM 2 M 3 . 

Let us now consider a (flexible) triatomic molecule. 
For simplicity, we freeze p\ and p 2 and consider the 
two-degree-of-freedom system consisting of pure bend- 
ing coupled with overall rotations. The potential V is 
then a function of cj> only and we choose it to be of the 
form shown in Fig. [3] This potential has two minima at 
4> = and cj> = it which correspond to two different lin- 
ear isomers which are separated by a barrier at <j> = it/2. 



This type of potential occurs, e.g., in the HCN/CNH iso- 
merization problem [10]. We consider the equilibrium 
which for a given magnitude J of the angular momen- 
tum arises from the barrier at <f> — tt/2 coupled with 
rotations about the first principal axis. In the absence 
of coupling between the bending and rotational degrees 
of freedom we would expect from the discussion of the 
rigid molecule above that this equilibrium is a saddle- 
center with a negative frequency w < 0. In fact also 
in the presence of coupling between the bending motion 
and the rotation this remains to be the case (at least for 
a moderate coupling strength). To illustrate the dynam- 
ical implication of this equilibrium we use the canonical 
coordinates (q,p) defined in ^ on the angular momen- 
tum sphere J 2 + J 2 + J 2 = J and construct a Poincare 
surface of section with section condition q = mod 2n, 
q > 0. Using the canonical pair ((j),p$) as coordinates on 
the surface of section we obtain Fig. |4j 

We see that, as expected, there appears to be a bar- 
rier associated with the momentum reaction coordinate 
Ptj,. Near = ir/2 no transitions are possible from the 
'reactants region' p^ > to the 'products region' p^ < 
for energies above the energy of the saddle whereas for 
energies below the energy of the saddle, transitions are 
possible. The rotational/vibrational motion of the tri- 
atomic molecule with the transition channel near 4> = tt/2 
being closed consists of unhindered bending motion be- 
tween the two isomers associated with cj> near and it, 
respectively, coupled with rotations. For motions with 
the channel near (j> = tt/2 being open, p^ can switch 
sign near (f> = tt/2 which corresponds to a trajectory 
bouncing back to the isomer it came from rather than 
switching to the other isomer. It is important to note 
that the saddle-center equilibrium studied here induces 
a local bottleneck for transitions between p^ > and 
p<f, < 0. Globally such transitions always occur in the 
present example when a trajectory passes close to the 
collinear configuration where <fi is close to or tt. This 
explains how a single trajectory can contribute points to 
the lower and the upper half in Fig.ga) even though the 
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FIG. 4: Poincare surface of section for the Hamiltonian |5| with 
parameters J = 0.2, p\ = 1 and p2 = 2 (see the text). Each picture 
is generated from a single trajectory of energy 0.0205625 in (a) and 
0.0197526 in (b), respectively. 



local channel near <f> — it/2 is closed. 

More generally: mixed positive and negative 
frequencies.- Near a saddle-center-- •• -center equi- 
librium the Hamiltonian can always be brought into the 
form ([!]). The NHIM at energy E is then obtained from 
the intersection of the center manifold of the equilibrium 
given by qo = po = and the energy surface H = E, i.e. 
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(pI 



h.o.t. = E- E 



(10) 



All studies on the geometric theory of reactions stud- 
ied so far [13] concern the case of positive frequencies 
Wfc. This corresponds to the Hamiltonian restricted to 
the center manifold having a minimum. If all frequen- 
cies are negative then the Hamiltonian restricted to the 
center manifold has a maximum. In this case the NHIM 
is again a (2/ — 3)-dimensional sphere. However as op- 
posed to the case of positive frequencies, po rather than 
q Q is the reaction coordinate as discussed in the present 
example for / = 2. Also the case of mixed signs of the Wfc 
occurs in many applications. It, e.g., also occurs in our 
example of the rotational-vibrational motion of triatomic 
molecule if we take the other vibrational degrees of free- 
dom into account which we for simplicity considered to 
be frozen. In the case of mixed signs of the Uk the NHIM 
is not a sphere but a non-compact manifold. Although 



no bottleneck in the energy surface is induced in these 
cases the NHIM has important dynamical implications 
as it has stable and unstable manifolds which are of one 
dimension less than the dimension of the energy surface 
and hence form impenetrable barriers. This bears some 
similarities to the case of non-compact NHIM's that have 
recently been considered for rank 2 and higher rank sad- 
dles [TU [15] . The case of mixed positive and negative 
frequencies will form an interesting direction for future 
research. 
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